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Abstract 

We study the model of interacting agents proposed by Chatterjee (2003) that al- 
lows agents to both save and exchange wealth. Closed equations for the wealth 
distribution are developed using a mean field approximation. 

We show that when all agents have the same fixed savings propensity, subject to 
certain well defined approximations defined in the text, these equations yield the 
conjecture proposed by Chatterjee (2003) for the form of the stationary agent wealth 
distribution. 

If the savings propensity for the equations is chosen according to some random 
distribution we show further that the wealth distribution for large values of wealth 
displays a Pareto like power law tail, ie P{w) ~ w^~^°'. However the value of a for 
the model is exactly 1. Exact numerical simulations for the model illustrate how, as 
the savings distribution function narrows to zero, the wealth distribution changes 
from a Pareto form to to an exponential function. Intermediate regions of wealth 
may be approximately described by a power law with a > 1. However the value 
never reaches values of ~ 1.6 — 1.7 that characterise empirical wealth data. This 
conclusion is not changed if three body agent exchange processes are allowed. We 
conclude that other mechanisms are required if the model is to agree with empirical 
wealth data. 
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1 Introduction 



The distribution of wealth or in come in society has b een of great interest 
for many years. Itahan economist Vilfredo Paretol (1897) was the first to sug- 
gest it followed a "natural law" where the higher end of the wealth distribu- 
tion is described by power law, P (w) ^ Repeated e r npirical studies 

Lew Solom(^ (|l997l ): lD^agulesciJ (|2nnih : lR,eed HughesI (|2nn2h : lAoyama Soumal 
show that the power law tail exhibits a remarkable spatial and tem- 
poral stability and while the value of the exponent, a, may vary slightly, it 
changes little from the value ~ 1.5. 



Even though the collected data stem from different sources and can be in- 
complete because of difficulties in accessibility (poor conclusions from income 
data in Sweden in iLevv SolomonI (|1997|) due to a too small number of wealth 
ranges in the data; total net capital of individual at death in the United States 
(US) reported to the Bureau of Cen sus and the Inland Revenue for tax her- 
itage purposes in iDragulescul (|2001[ ): distributions of sizes of incomes, cities, 
internet files, biological taxa, gene family and protein famil y frequencies in 
R.eed Hughes! (|200a); and income distributions in the Japan in lAovama Soumal 
the common conclusion which can be drawn is that the high end that 
exhibits the power law is characterised by several multiples or even tens of 
multiples of the average income/wealth (only 5% of population income-data 
in the US conforms to a power-law and the power la w for the yearly inc ome 
data in the United Kingdom sets in only for > 50k £ iDragulescd (|2001l ). in- 
come distributions in the Japan in 2000 exhibit power laws only for > 5 • 10^ 
thousands of Yen) . 



For around 100 years the tantalising Pareto law remained without explana- 
tion. The renewed interest by physicists and mathematicians in econo- and 
sociophysics has however led to publicat ion of a number of new papers on the 



topic in recent years (see ISlaninal (|2004f l for an extensive literature review). 



The fact that multiplicative power law processes can lead to power law distri- 
butions has beei i known for many years from st udies as divers e as the frequency 
of words in text lYuld (|l997l ). eco nomi c grow th iGibrIi] ^93^, city populations 
Zip! (Il949f). w ealth distribution liiri (1977) and stochastic renewal processes 
KestenI f)l973h . 



In the analysis of these distributions ISolomon ( 2002[ ) has recently proposed the 
use of Generalised Lotka Volterra (GLV) equation that combines a multiplica- 
tive random process with an autocatalytic process. The latter redistributes 
a fraction of the total money to ensure the money possessed by an agent is 
never zero. This simulates in a simplistic way the effect of a tax. The model 
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equations lead to a wealth distribution P{w) of the form: 

g(l-o)/«) 



P w) 



where and a — 1 is a positive number that is a ratio of parameters of the model 
that are related to social security and some random investements respectively. 
For large values of income w this indeed exhibits a Pareto behaviour. 

However two issues arise. The first is that empirical studies of income dis- 
tributions show that this function does not describe wel l the very low en d 
of the income distribution which is essentially exponentia JPragulescu ( 2001 ). 



The second relates to use of the multiplicative stochastic term. It is certainly 
necessary to secure the right form for the distribution function but how does 
it arise in the first place? 



More recentlv IChatteri have developed a model of pairwise interact- 

ing agents i and j that exchange money by analogy with an ensemble of gas 
molecules that exchange momentum. In Chatteriea ( 2003[ l's model, however, 



the agents are allowed to save a fraction Aj of their money prior to an in- 
teraction. The total money held between two agents is conserved during the 
interaction process. The governing equations for the evolution of wealth Wi 
and Wj of agents i and j respectively are given by: 

Wi (t + l) = XiWi + e [(1 - Xi)wi + (1 - \j)wj] 

Wj (t + 1) = XjWj + [1 - e] [(1 - \,i)wi + (1 - \j)wj] (2) 



Here each agent, i, has a savings propensity, Aj. The remaining money is 
divided during the exchange process in a random manner determined by a 
uniformly distributed random nurnber e between zero and one. From their 
numerical calculations Chatter ie3 (|2003[ ) found the following results for the 



stationary wealth distribution P{w). Here 

(1) With no saving (Aj = for all i) agents behave randomly and the dis- 
tribution follows the Gibbs rule P{w) ~ exp{—w/ (w)) where (w) is the 
average wealth of agent. P{w) has a maximum when w = 0. 

(2) If the saving propensity is non-zero and takes the same constant value 
for all agents (Aj = A for all i) the resulting distribution can be fitted 
well by the heuristic function: 

71 ^ 

P (w) = ^^w""-^ exp i-nw) (3) 
Tin) 

where r(n) is the gamma function and the parameter n is related to the 
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saving propensity, A as follows: 



The power, w"'~^, qualitatively changes the distribution so that it has a 
maximum for w > 0. The author does not give any theoretical arguments 
for the use of this distribution. 
(3) If the saving propensity for the agents is chosen according to some ran- 
dom distribution, like uniform or power-law distributions with < A < 1, 
the numerical output for large values of money gives P (w) ~ w^^~°'. 
This is the celebrated Pareto law. Numerical calculations yield a value 
for a = 1.03 ± 0.03. The authors show empirical data for wealth distri- 
butions for both Japan and the USA. These data clearly exhibit power 
laws with values of a greater than this value. However the authors leave 
the reader wondering whether the model could fit this data better. A 
further calculation allowing only a fraction, p, of agents to save is made 
but the the value for the Pareto law remains unchanged. The authors 
do not investigate the possible changes in a as a result of using savings 
distributions that differ from uniform distributions. 



This work is interesting in that it brings together within one framework the 
distributions of both Gibbs and Pareto. However it leaves open tantalizing 
questions. 

(1) Is it possible to predict analytically the expressions 3 and 4? 

(2) How does the value of a within the model depend on the nature of the 
savings distribution? 

(3) Could it be that the value of a is actually unity? 

(4) Is there a way of reconciling the approach based on the GLV equations 
and the exchange theory of Chakrabarti? 

Note that caution has to be taken by fitting the models described above to 
incom e data obtained fro m Inland Revenues in different countries. As pointed 
out in iDragulesciJ (|200d ) the wealth has to be understood as a commodity 



that is subject to an incessant process of exchange rather then as valueables 
like precious metals, "hard currency", bonds or works of art that have been 
deposited in a bank account in order to serve as a lifetime security. In this sense 
the distributions that come out of the models should be fitted to a momentary 
distribution of money in the society; a distribution they may or may not 
be in equillibrium. Since people rarely disclose their momentary wealth the 
statistical data one avails of regards more the total wealth of individuals, ie 
the wealth that has been accumulated throughout their whole lifetimes and 
is reported to the Revenue office only after death (to fulfill the heritage tax 
requirements). Since, however, individuals with small and medium wealths are 
rather unlikely to invest parts of their income in any sort of lifetime securities. 
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because their earnings are small and are spent in their total to cover the cost 
of living, the low end of the momentary money distribution in equillibrium 
should coincide with the low end of the wealth distribution obtained from the 
Revenue data. Differences will only be observed in the high end. 



In the next section we develop the theory for the model by Chatter ie3 ( 20031 ) 



and show that it is indeed possible to demonstrate that the conjecture sum- 
marised above is, to within a certain well defined approximation, correct. We 
then study asymptotically the behaviour of the wealth distribution for the 
case where the savings propensity varies for the agents and demonstrate that 
if the wealth distribution P{w) ~ w^~^°' then a is exactly unity for this model. 
We demonstrate in section 2.2 that this conclusion remains unchanged even 
when three agent exchange processes are allowed. The conclusions of the mean 
field analytic analysis are supported by exact numerical simulations shown in 
Figs. 5 and 6. 



2 Theoretical analysis 



Complete information about the processes at time t is given by the N agent dis- 
tribution function fN{vi, . . . ,vn)- In what follows we shall assume the mean- 
field approximation. This implies that the A^-agent distribution function: 

Mvi, ...,vn)=P (^f]ik <Vi<Vi + dv^ /{dvf = Y[fi{vi) (5) 



We can now invoke the Boltzmann equation lErnstI (|l98l[ ) for the one-agent 
wealth V distribution fi{v;t) at time t. Thus: 



dji{v; t) = / dwdv'dw' {W {v w\v' w') h{v')fi{w') - W {v' w'\v w) h{v)fi{w)) (6) 



where the transition probabilities W () are given via the rules for the collision- 
dynamics (2): 



W {v' w'\v w) = 

6 {v' - {Xv + e(l -X){v + w))) • 5 K - {Xw + (1 - e)(l - X){v + w))l7) 

Introducing the Laplace transform fi{x;t) := fi{v]t) exp{—vx)dv we ob- 
tain an integro-differential equation for the temporal evolution: 

dtfiix; t) + ~h{x- 1) = (/2 (Ax + e (1 - A) X, e (1 - A) x; t)) (8) 
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where the random spatial variability of the saving propensities A and exchange 
fractions e is accounted for by the averaging process () over their random 
distributions. 



We note at this point that this model assumes elastic scattering, i.e. conserva- 
tion of wealth during the exchange process (2) and the existence of a stationary 
solution. This is in contrast to many previo us models f ormul a ted in different 
contex t s where 'wealth' may be eith e r lost IKrapivskv l (EooS); Ib en-AvrahamI 



( 20031) : iBobvlev (2000,); iBaldassarril (|2002l) or gained Slanina (200 j) in the 



exchange process and the distribution function has a power-law-tail only in 
an asymptotic sense. 

We now write the stationary solutions = lirrit >oofi{x]t) both in terms 

of solutions of non-linear integral equations (Master Equations (MEs)) and 
in terms of expansions fi{x) = Y^'^=o{—^)""mnX^ over moments (Ms) (f") = 
rUn ■ n\ which satisfy recursion relations. It is convenient to distinguish two 
cases: 

(I) Saving propensity: A 7^ but equal for all agents: 



(l-A)x 



ME: xfi{x) 



1 - A 



f\{Xx + 4>)f\{(P)d(P 



(9) 



Ms: TTLp = ^ mgmp_g Ci''\X) with 
and c'a^^'""-''^^'^"" 



p-q 



1-A 

^ with C(^) = ^i^ 
P+l 



(II) Random saving propensity: A ~ Pa (A). 

ME: xh{x) = jd\P^ I A(Aa; + 0)/i(0)ci0 




Here we waive the writing of equations for the moments since due to to the 
wealth distribution having a power-law tail they may not of course exist. 
Now an assumption about an asymptotic expansion in the "wealth-domain" : 
fii'^) = Y.'^=o(^n/v^~^°'~^^ leads to a decomposition of the function in the 
Laplace domain into two parts 

fiix)=fr^ix)+fr^ix) (12) 
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with the first part being an analytic function fi^^{x) = 1 — x + O (a;^) and 
the second part fi^^^{x) = Y.'^=o bnx"''^" having a leading term x" of order a. 

2.1 Conjecture by Patriarca, Chakraborti and Kaski (PCK): 



Solving the moments' equations (10) with initial conditions mo = 1 and mi = 
1 recursively, ie. expressing, via the pth equation, function of A and 

all previous values of m, (ie mo, mi, ... , mp_i), one obtains: 



A+2 A+2 

^^=2(rT^ "^^ = 2(rTW ^ ^ 

72 + 12A - 2A2 + 9A3 - 
~ 24(1 + 2A)2(3 + 6A - A2 + 2A3) ^ ' 

The first three moments mi,m2 and ms coincide with the moments of PCKs 
function (3) if the relation between the parameters n and A is given by (4). 
Indeed the coefficients of a series expansion of the Laplace transform 

oo 

n 



P{x) = / P(0 exp(-ex)rfe = —— (15) 

\x + n/ 







of the function (3) agree with moments (14) up to the third order subject 
to equation (4) being satisfied. This is shown in a nice way in Fig. 1. The 
deviation A/i(a;) between the exact solution of the ME (9) and the ansatz (3) 
has a leading fourth order: 

' 8n3(10n3 + 30n2 + 45n-4) + ^ l"^ J 



It is hard to say if a more general class of functions than (3) would satisfy the 
ME to higher expansion orders. 



2.2 The power-law tail: 



Calculations by Chatter ie3 ( 20031 ) and ourselves suggest that the value of the 



exponent a is equal to one. Let us look at this aspect in more detail. Inserting 
the expansion (12) into the ME (9) and comparing coefficients of order x°'~^^ on 
both sides of the equation leads to a transcendental equation for the exponent, 
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A{X; n) 
0.25 



0.2 



0.15 



0.1 



0. 05 



Num. n = 2 

Num. n = 3 

Num. n = 4 

Num. n = 5 




1/4 



0.4 0.6 
2/5 3/6 4/7 



Fig. 1. Deviations A(A;n) = J2l=o ("^pC-^) — rrip""-^ (n)) / p\ of the exact moments 

rUp of the wealth distribution from the moments rrip^^^ derived from the conjecture 
plotted as a function of A for n = 2, 3 . . . , 9. Solid lines (dot symbols) correspond to 
analytical (numerical) solutions of the moment equations (10). We see that the min- 
ima A = (n- l)/(n+2) = {1/4,2/5,3/6,4/7,5/8,6/9,7/10,8/11} of the deviations 
do correspond to the PCK conjecture (4). 



a: 



((1 - A)") + 



1 - A 



a+l ' 



A 



= q; + 1 



(17) 



Clearly if we choose a = 1, (17) we obtain an indentity for any distribution 
of A. This would seem to be true even for a distribution that assumes only 
a fraction p of the agents save and the remainder do not save, i.e. Pa ''(A) = 
PPa(A) + (1-p)<5(A). 

However, whether other solutions for a exist is an open question and depends 
on the distribution of the saving propensity Pa (A). We try to clarify this ques- 
tion below. 

For uniformly distributed propensities Pa(A) = I//2 for < A < /2 < 1 the 
only solution is a = 1 (see Fig. 2). Likewise if Pa (A) is a normal distribution 
with a variable mean I2 where I/2I < 1 and the standard deviation is small 
(See Fig. 3) or with a fixed mean and variable standard deviation (Fig. 4) the 
exponent similarly turns out to be unity. 

Now we make a stronger statement and say that there is no continuous and 
differentiable distribution of saving propensities A G [0, 1] that would yield 
a ^ 1. Indeed since every distribution Pa(A) can be constructed as a weighted 
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Fig. 2. The left-hand-side of equation (17) plotted as a function of a for a uniformly 
random saving propensity A with < A < /2 for different values of I2 (dashed lines) 
and the right-hand-side (rhs) of the equation (17) (solid line). As we can see there 
is no other solution of the transcendental equation (17) except a = 1 in the range 
a € [0,2]. 

(possibly continuous) linear combination of uniform distributions 
1 

p^{X) = J w{u)U{0,u)du (18) 




and since for a uniform distribution U (0, u) the left-hand side of the transcen- 
dental equation (17) intersects the right-hand side only for a = 1 in the range 
a G [0,2] then the last statement holds also for a generic distribution Pa(A). 
Here we used conditional averaging. This means that the average in equation 
(17) is carried out as an average over U{0, v) conditioned on v first and then 
over the distribution w{v) of random values of v. 



2.3 Beyond the mean- field approximation: 



Many-agent distribution functions fN^^i, . . . ,X]\f) may not be produced cor- 
rec tly within the mea n-field approach. Furthermore the wealth-exchange model 
by Chatter ie3 ( 2003[ l may be extended to N-point interactions: 



(t + 1) = XiWi + ti 



N 

E( 



(19) 
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Fig. 3. The same as in Fig. 2 but for the propensity A conforming to a truncated 
|A| < 1 normal distribution with variable mean I2 and standard deviation 0.01. 
As before dashed lines denote the left-hand side and the solid line denotes the 
right-hand side of equation (17). Again the only solution of the transcendental 
equation is a = 1 in the range a € [0, 2] . 




Fig. 4. The same as in Figs.2 and 3 except that now the mean of the normal 

distribution of propensities A is fixed and equal to 0.5 and the standard deviation 
varies. Here again no new solutions except a = 1 are obtained. 

Here we mean that at every time step exchange processes involving any number 
of agents can happen - each with a certain likelihood. We perform the analysis 
ior N — 3 in order to find out what kind of mathematical difficulties we will 
come across. Now the master equation for the 2- agent distribution function in 
the Laplace domain (compare with (8)) reads: 
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dtf2(x,y;t) + f2(x,y;t) = 

(1 - a) (/2 (Xx + Ai (ex + (1 - e)y) ,Xy + Xi (ex + (1 - e)y) ; t)) + 

^ ^/a (Ax + Ai (eix + e2y) ,Xy + Xi (eix + esi/) , Ai (tix + est/) ; t)) + 
^ (/s (Ai (€23; + esy) , Ax + Ai (e2a; + e^y) ,Xy + Xi (e2X + e^y) ; tfj (20) 

where A + Ai = l,e + ei<l and a and (1 — cr) denote likelihoods of three- 
agent and two-agent exchange processes respectively. The first (second and 
third) term(s) on the right-hand side in (20) account(s) for two- (three-) agent 
exchange processes repectively. Setting |/ = wc obtain equation (8) except 
for three-agent exchange terms that were neglected in the first place and now 
have been added appropriately. Setting a; = y = we obtain an identity from 
the normalisation condition /2(0,0) = /3(0,0,0) — 1. 



2.4 The power-law tail with three-agent exchange processes: 



Now the transcendental equation, derived from the master equation (20), has 
the following form: 



[I -a) 
2 



q; + 2 



((1 - A)") + ^ 
2((1-A)") + 



1 - A 



a+l ' 



+ 



1- A 

/l-(a + 2)A"+^ + (a + l)A"+^' 
\ (1 - A)^ 



= a + 1 (21) 



and a q; = 1 is again the only solution (compare Fig. 5) of this equation for ar- 
bitrary saving propensity distributions Pa (A) and for any likelihood a G [0, 1] 
of three-agent exchange processes. This is in conformance with our numeri- 
cal simulations that also show that an introduction of three-agent exchange 
processes do not alter the exponent. 



2.5 Moment equations in the case of two- and three-agent exchange processes: 



The expansion of the steady-state solution in terms of two-agent correlations 
{v^w'^) = rrip^q ■ (p\q\) now reads f2(x,y) = Zl^g=o(~l)^^'^^p,5^^2/^ the 
two- agent correlations satisfy following equations: 
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Fig. 5. The second term (corresponding to three agent exchange processes) on the 
left-hand-side of equation (21) plotted as a function of a for a uniformly random 
saving propensity A with < A < ^2 for different values of Z2 (dashed lines) and 
the right-hand-side (rhs) of the equation (17) (solid line). We see that three-agent 
exchange processes do not lead to a change of the exponent a. from unity to a 
different value. 



m. 



p,<i 




"^pi +P2,qi +92 ^^^^ 



n^pi +P2 ,qi+q2,ri +r2 -^pi 



PI + 11 + ri 
V2 + 12 + r2 = q 



where 



(^^:£(A) = ^ / (A + ^-'^ (1 - (Ai - o'-'^ 




(23) 



If we neglect ternary exchange processes {a — 0) we end up with following 
equations for the moments of the 1- agent distribution function: 



1712 = rn^i 



A + 2 
2(1 + 2A) 



1 + 2A 



12 



(4 - 2A + 2A2 + A3) -6 + 3A + 2X^ + A^ 

- '"^'^2(3 + 6A-A2 + 2A3) ~ ""^'^6(3 + 6A - A^ + 2A3) 

_ ((2 - 2A + 3A^)mi,4 - 2(-l + A)m2,3) 

(4 + 8A-3A2 + 6A3) ^^^^ 

These equation reduce to equations (14) under the mean- field approximation 

m„n = m„ma. 



3 Exact computer simulations of the model 



Intuitively, it is clear that the distribution function depends critically on the 
relative values of the mean value of the savings propensity and the spread or 
mean square deviation of the saving distribution function. As the spread of 
saving propensities tends to zero, the distribution function must change from 
one having a power law tail to one with an exponential tail. How does this 
change take place? Could it be that the effective power law region shifts to 
take on values of a > 1? Might this explain the empirically observed facts? 
We have made computer simulations of the model for different values of these 
parameters. The results are described below. 



A uni form distribution of the savings parameter A in the model of Chatterie3 



results in a power law distribution of the cumulative distribution shown 



in the Fig. 6. 

The situation is different in the case of A being Gaussian distributed. Here 
the cumulative wealth distribution may still be approximately described by a 
power law for widths between 1 and 0.45 (see Fig. 7). 

One thing seems clear. The region where power law behaviour is observed is 
fairly well marked even when the spread is small. And moreover the slope is 
consistent with higher values of a. However it does not seem that a can take 
on values greater than around 1.2. As it stands then we conclude that this 
model does not offer a complete picture of the empirically observed wealth 
dynamics. 



4 Conclusions 



We have studied the model of interacting agents proposed bv lChatteri e3 (|2003l) 
that allows agents to both save and exchange wealth. Closed equations for the 
wealth distribution are developed using a mean field approximation. 
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Fig. 6. Results of computer simulations, obtained for 500 agents, and 39000 reali- 
sations, each taken after 100000 equilibration steps. The data is well described by 
a power law with exponent -1 in the range of 0.3 up to about 15 times the average 
money per agent. Changing the width of the distribution of A immediately leads 
to a loss of the power law. The latter may be a finite size effect since our analyt- 
ical calculations from section 2.2 show that the power-law does occur for uniform 
distributions of any widths. 



We have shown that when all agents have the same fixed savings propensity, 
subject to certain well defined approxin iations defined in the text, these equa- 
tions yield the conjecture proposed bv IChatteri e3 (Hooi) for the form of the 
stationary agent wealth distribution. 



If the savings propensity for the equations is chosen according to some random 
distribution we have further shown that the wealth distribution for large val- 
ues of wealth displays a Pareto like power law tail, ie P{w) ~ w^~^°'. However 
the value of a for the model is exactly 1. Exact numerical simulations for the 
model illustrate how, as the savings distribution function narrows to zero, the 
wealth distribution changes from a Pareto form to to an exponential function. 
Intermediate regions of wealth may be approximately described by a power 
law with a > 1. However the value never reaches values of ~ 1.6 — 1.7 that char- 
acterise empirical wealth data. This conclusion is not changed if three body 
agent exchange processes are allowed. We conclude that other mechanisms are 
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Fig. 7. The same as in Fig. 6 but for the saving propensity being Gaussian dis- 
tributed. Here the power law exponent decreases with decreasing width, but only 
over a small range from -1 to -1.13. Narrower Gaussian distributions do not re- 
sult in wealth distributions that can be described by power law. (Note that in the 
limit of small widths both Gaussian and uniform distribution give the same wealth 
distribution.) 

required if the model is to agree with empirical wealth data. 
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